Controllability of ferromagnetism in graphene 
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We systematically study magnetic correlations in graphene within Hubbard model on a honey- 
comb lattice by using quantum Monte Carlo simulations. In the filling region below the Van Hove 
singularity, the system shows a short-range ferromagnetic correlation, which is slightly strengthened 
by the on-site Coulomb interaction and markedly by the next-nearest-neighbor hopping integral. 
The ferromagnetic properties depend on the electron filling strongly, which may be manipulated by 
the electric gate. Due to its resultant controllability of ferromagnetism, graphene-based samples 
may facilitate the development of many applications. 



The search for high temperature ferromagnetic semi- 
conductors, which combine the properties of ferromag- 
netism (FM) and semiconductors and allow for practi- 
cal applications of spintronics, has evolved into a broad 
field of materials sciencep], Q. Scientists require a ma- 
terial in which the generation, injection, and detec- 
tion of spin-polarized electrons is accomplished with- 
out strong magnetic fields, with processes effective at or 
above room temperature [3|. Although some of these re- 
quirements have been successfully demonstrated, most 
semiconductor-based spintronics devices are still at the 
proposal stage since useful ferromagnetic semiconductors 
have yet to be developed 0. Recently, scientists antic- 
ipate that graphene-based electronics may supplement 
silicon-based technology, which is nearing its limits 0, Q. 
Unlike silicon, the single layer graphene is a zero-gap 
two-dimensional (2D) semiconductor, and the bilaycr 
graphene provides the first semiconductor with a gap 
that can be tuned externally 0- Graphene exhibits gate- 
voltage controlled carrier conduction[7l-[loj . high field- 
effect mobility, and a small spin-orbit coupling, making 
it a very promising candidate for spintronics application 
[ill , [lil ]. In view of these characteristics, the study of 
the controllability of FM in graphene-based samples is of 
fundamental and technological importance, since it in- 
creases the possibility of using graphene in spintronics 
and other applications. 

On the other hand, the existence of FM in graphene 
is an unresolved issue. Recent experimental and theo- 
retical results in graphene[l~3- 15 1 show that the electron- 
electron interactions must be taken into account in order 
to obtain a fully consistent picture of graphene. The 
honeycomb structure of graphene exhibits Van Hove sin- 
gularity (VHS) in the density of states (DOS), which may 
result in strong ferromagnetic fluctuations, as demon- 
strated by recent quantum Monte Carlo simulations of 



the Hubbard model on the square and triangular lat- 
tices [l6, 17 1 . After taking both electron-electron inter- 
action and lattice structure into consideration, the bidi- 
mensional Hubbard model on the honeycomb lattice [3- 
l2fjj is a good candidate to study magnetic behaviors in 
graphene. Early studies of the bidimensional Hubbard 
model on the honeycomb lattice were based on mean field 
approximations and the perturbation theory [2(|. How- 
ever, the results obtained are still actively debated be- 
cause they are very sensitive to the approximation used. 
Therefore, we use the determinant quantum Monte Carlo 
(DQMC) simulation technique 2l|, l22| to investigate the 
nature of magnetic correlation in the presence of moder- 
ate Coulomb interactions. We are particularly interesting 
on ferromagnetic fluctuations as functions of the electron 
filling, because the application of local gate techniques 
enables us to modulate electron filling [7H10(, which is the 
first step on the road towards graphene-based electronics. 

The structure of graphene can be described in terms of 
two interpenetrating triangular sublattices, A and B, and 
its low energy magnetic properties can be well described 
by the Hubbard model on a honeycomb lattice 18 - 2(| , 
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where t and t' are the nearest and next-nearest-neighbor 
(NNN) hopping integrals respectively, fi is the chemical 
potential, and U is the Hubbard interaction. Here, a ia 
( a la) annihilates (creates) electrons at site Ri with spin 
u (ct=|,|) on sublattice A, b ia (b\ a ) annihilates (creates) 
electrons at the site Ri with spin a (c=t, -i) on sublattice 
B, n aia =a\ a a icr and n Ua =b\ a b ia . 

Our main numerical calculations were performed on 
a double-48 sites lattice, as sketched in Fig. [TJ where 
blue circles and yellow circles indicates A and B sub- 
lattices, respectively. The structure of the honeycomb 
lattice leads to the well known masslcss-Dirac-fcrmion- 
like low energy excitations and the two VHS in the DOS 
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<n>=0.25 



FIG. 1: (Color online) (a) Sketch of graphene with double-48 
sites; (b) First Brillouin zone and the high symmetry direction 
(red line); (c) DOS (solid dark lines ) and fillings (n) (dash red 
lines ) as functions of energy with t'=0.10t; and (d) t'=0.20t. 



(marked in Fig. 1) at <?i>=0.75 and 1.25 correspond- 
ing to E=—2t' ± t, respectively as t' < t/6. While when 
t' > t/6, a third VHS appears at the lower band edge, 
which is a square root singularity marking the flattening 
of the energy band near T point. They determine much 
of system's properties. According to the values of t and 
U reported in the literatures for graphene[2^ 
the ratio U/\t\ maybe expected to be 2.2 ~ 6.0 
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which 

is around the range of half- bandwidth to band width pjj. 
where the mean filed theory does not work well while the 
DQMC simulation is a useful tool [22[. The exact value 
of t' is not known but an ab initio calculation [26} found 
that t' It ranges from 0.02 to 0.2 depending on the tight- 
binding parameterizations. Therefore, it is necessary to 
study the ferromagnetic fluctuations within the Hubbard 
model on the honeycomb lattice by including t 1 . 

In the followings, we show that the behaviors of mag- 
netic correlation are qualitatively different in two filling 
regions separated by the VHS at <n>=0.75. In the filling 
region below the VHS the system shows a short-ranged 
ferromagnetic correlation and the on-site Coulomb in- 
teraction tends to strengthen ferromagnetic fluctuation. 
The ferromagnetic properties depend on the electron fill- 
ing, which may be manipulated by the electric gate. Fur- 
thermore, the ferromagnetic fluctuation is strengthened 
markedly as t' increases. Our results highlight the cru- 
cial importance of electron filling and the NNN hopping 
in graphene. The resultant controllability of FM may fa- 
cilitate the new development of spintronics and quantum 
modulation. 

To study ferromagnetic fluctuations, we define the spin 
susceptibility in the z direction at zero frequency, 



X(q) = dr 
Jo 
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where m !o (T)~e HT m,i a (0)e Ht with m !a =aLaff — o/^a^ 
and m,i b =bl^bit — b\^bi^. Here x is measured in unit of 




FIG. 2: (Color online) At <n>=0.25, inverse of magnetic 
susceptibility, l/x(5=r) versus temperature with (7=3|t|, 
t'=t/W, t/6, and t/5. Fitted line l/x(r)=a(T - 6) are also 
shown. Inset: Magnetic susceptibility x(l) versus q at differ- 
ent temperatures with t'—0.10t and t/=3|t|. 




FIG. 3: (Color online) Magnetic susceptibility x(l) versus 
momentum q at different t , here U—3 |t|, <n>=0.25 and 
T=0.10|t|. 



\t\ , and x(r) measures ferromagnetic correlation while 
x(K) measures antifcrromagnetic correlation. 

We first present temperature dependence of the mag- 
netic correlations at <n>=0.25 with different t' and U. 
Fig. [2] shows l/x(g=r) versus temperature at U=3|t| 
with t'=t/10, t/6, and t/5. Data for U=5\t\ as t'=t/5 
are also shown. In the inset, we present x{q) versus mo- 
mentum q at different temperatures with t'—t/lQ and 
C/=3|t|. It is obvious that xil) A as strong temperature 
dependence and one observes that x(M) and x{K) grow 
much slower than x(r) with decreasing temperatures. 
Moreover, l/x(r) exhibits Curie-like behavior as tem- 
perature decreases from \t\ to about 0.1 \t\. Fitting the 
data as l/x(F)=a(T - 9) (solid lines in Fig|2]) shows 
that 9 is about 0.02|t|~ 580^ at t'=t/5, and we also 
note that both O and x(r) are enhanced slightly as the 
on-site Coulomb interaction is increased. Positive values 
of indicate that the curves of l/x(r) start to bend 
at some low temperatures and probably converge to zero 
as T — > 0, i.e., x(T) diverges. This demonstrates the 
existence of ferromagnetic state in graphene. 

From Fig. [5J we may also notice that t' plays a re- 
markable effect on the behavior of x(o)i an(1 results for 



3 



C 4 
ii 

15 2 



(a) 

U=3|t| -°- x(0 
T=0.167|t| X (K) 

f=o.iot — x*ryx(K) 


(b) /\ 

f=0.20t P \ 









0.0 0.2 0.4 0.6 0.8 0.0 0.2 0.4 0.6 0.8 1.0 

<n> <n> 



FIG. 4: (Color online) Magnetic susceptibility x( a= ^) (red) 
and x(°~ K) (dark) versus electron filling at C/=3|t| and 
T=0.167|t| with (a) t'=Q.lt and (b) t'=Q.2t. 



X{q) dependent on q with different t' at £7=3|i|, T=0.10|i| 
and <n>~0.25 have been shown in Fig. [3] Clearly, x(r) 
gets enhanced greatly as t' increases, while xiM) an d 
x{K) increase only slightly. Thus, again it is significant 
to demonstrate that ferromagnetic fluctuation gets en- 
hanced markedly as t' increases. Furthermore, the strong 
dependence of FM on t' suggests controllability of FM in 
graphene by tuning t' . 

A great deal of current activity in graphene arises from 
its technological significance as a semiconductor material 
where carrier density can be controlled by an external 
gate voltage 0-0] • To understand filling dependence of 
magnetic correlations intuitively, we present x(F) (red), 
x{K) (dark), and their ratio x(T)/x(K) (blue) versus 
filling for (a) t'=0.lt and (b) t'=0.2t in Fig. gj where 
U=3\t\ and T=0.167|i|. There is a crossover between 
x(r) and x(K)i which indicates that the behaviors of 
x(q) ar c qualitatively different in two filling regions sep- 
arated by the VHS at <n>=0.75, where the ratio is 1. 
This is due to the competition between ferromagnetic 
and antiferromagnetic fluctuations. The antiferromag- 
netic correlations are strong around the hall-filing and 
they may dominate the shape of x(q) m a wide filling 
range up to the VHS. The effect of t' in enhancing ferro- 
magnetic fluctuation can also be seen by comparing Fig. 
|4](a) with (b). At electron filling <?i>=0.25, the ratio of 
x(T)/x(K) is about maximum for t'=0.2t and is substan- 
tial for i'=0.1i, which is the reason why did we choose 
electron filling 0.25 in Figs. [3J and EH From the global 
picture shown in Fig. g] it is clear that the strength of 
ferromagnetic correlation strongly depends on the elec- 
tron filling, which may be manipulated by the electric 
gates in graphene, since n oc V^[7H9|- The filling region 
for inducing FM required here likely exceeds the current 
experimental ability. In fact, the challenge of increasing 
the carrier concentration in graphene is indeed very im- 
portant and it is a topic now in progress. The second 
gate (from the top) and/or chemical doping methods are 
devoted to achieving higher carrier density [l(J. More- 
over, our results present here indicate the electron filling 
markedly affects the magnetic properties of graphene, 



and the controllability of FM may be realized in fer- 
romagnetic graphene-based samples. Furthermore, the 
change of ferromagnetic correlation with t' may also lead 
to controllability of FM in graphene. For example, one 
can tune t' by varying the spacing between lattice sites. 
Tuning t' can also be realized in two sub triangular op- 
tical lattices 12 "i 
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Due to the peculiar structure of 
graphene, which can be described in terms of two inter- 
penetrating triangular sublatticcs, making controlling t' 
in principle possible in ultracold atoms system by using 
three beams of separate laser [29j. 

In summary, we have presented exact numerical results 
on the magnetic correlation in the Hubbard model on a 
honeycomb lattice. At temperatures where the DQMC 
were performed, we found ferromagnetic fluctuation dom- 
inates in the low electron filling region, and it is slightly 
strengthened as interaction U increases. The ferromag- 
netic correlation showed strong dependence on the elec- 
tron filling and the NNN hoping integral. This provides a 
route to manipulate FM in ferromagnetic graphene-based 
samples by the electric gate or varying lattice parame- 
ters. The resultant controllability of FM in ferromagnetic 
graphene-based samples may facilitate the development 
of many applications. 
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